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Abstract 

Using quasipotential approach, we have studied exclusive semileptonic decays of heavy mesons 
with the account of relativistic effects. Due to more complete relativistic description of the s quark 
more precise expressions for semileptonic form factors are obtained. Various differential distributions 
in exclusive semileptonic decays of heavy mesons are calculated. It is argued that consistent account 
of relativistic effects and HQET motivated choice of the parameters of quark-antiquark potential allow 
to get reliable value for the ratio A2 (0)/Ai(0) in the D — > K*lv\ decay as well as the ratio Y(D — > 
K*lvi)/T(D — > Klv{). All calculated branching ratios are in accord with available experimental data. 
q ; PACS number(s): 12.39K1, 13.20He, 13.40Hq 
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1 Introduction 



Semileptonic decays of heavy mesons provide an important tool to investigate quark dynamics and to 
Qh! determine Cabibbo-Kobayshi-Maskawa (CKM) matrix elements. Hadron dynamics is contained in form 

factors, which are Lorentz invariant functions of g 2 , the square of momentum transfer. These form factors 
cannot be calculated from the first principles of QCD by now. Thus various potential models, sum rules 
and lattice calculations have been proposed Recently considerable progress has been achieved in 

describing heavy meson decays by the use of heavy quark effective theory (HQET) |?J . It has been found 
that in the limit of infinitely heavy b and c quarks their mass and spin decouple from the dynamics of 
the decay and the description of a process such as B — ► Dlvi is strongly simplified. For D decays HQET 
predictions are less useful, because in this case symmetry breaking corrections appear to be rather large. 
It is also important to note that since B and D mesons contain light quark, relativistic effects are quite 
significant and consistent relativistic description of heavydight quark system is necessary. 

Our relativistic quark model (RQM), has some features that make it attractive and reliable for the 
description of heavy mesons. Firstly, RQM provides a consistent scheme for calculation of all relativistic 
corrections and allows for the heavy quark 1/toq expansion. Secondly, it has been found || that the 
general structure of leading, next- to- leading and second order 1/rag corrections in RQM is in accord with 
the predictions of HQET. The heavy quark symmetry and QCD impose some rigid constraints on the 
parameters of the long-range confining potential of our model. It gives an additional motivation for the 
choice of the main parameters of RQM and leads us to the conclusion that the confining quark-antiquark 
potential in meson is predominantly Lorentz- vector (with the Pauli term), while the scalar potential is 
anticonfining and helps to reproduce the initial nonrelativistic potential. This model has been applied to 
the calculations of meson mass spectra || , radiative decay widths |1(J , pseudoscalar decay constants [jllj , 
rare radiative |Q and nonleptonic [^3| decay rates. Semileptonic decays of B and D mesons have been 
considered in our model in fll4j . Here we refine our previous analysis with more complete account of 
relativistic effects and HQET constraints. We also consider exclusive decay spectra and q 2 dependence 
of form factors. 
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In Sect. 2 we briefly describe RQM. Sect. 3 is devoted to the calculation of form factors and semileptonic 
branching ratios and contains analytical expressions and numerical results for the differential distributions 
for the decays into pseudoscalar as well as vector final states. We give our conclusions in Sect. 4. 



2 Relativistic Quark Model 



Our model is based on the quasipotential approach in quantum field theory 15 . A quark-antiquark 
bound system with the mass M and relativistic momentum p in the center of mass system is described 
by a single-time quasipotential wave function ^Sm(p)i projected onto positive-energy states. This wave 
function satisfies the quasipotential equation 



M 2 - (p 2 



,2U/2 



-(P 2 



,?W 2 



*m(p) = 



d 3 q 

(2tt) 3 



K(p,q;M)* M (q), 



The quasipotential equation pi) can be transformed into a local Schrodinger-like equation p6| 



b 2 (M) 



*m(p) 



d 3 q 

(2tt) 3 



y(p,q;M)* M (q), 



where the relativistic reduced mass is 



E\E 2 
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4M 3 
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Ex = 



M 2 - mi 



E 2 = 



M 2 - m\ 



2M 2M 
and the square of the relative momentum on the mass shell is 



E 1 +E 2 = M. 



b 2 {M) 



[M 2 - (mi + m 2 ) 2 } [M 2 - [m x - m 2 ) 2 ] 



AM 2 



(4) 



mi )2 are the quark masses. 

Now it is necessary to construct the quasipotential V(p, q; M) of the quark-antiquark interaction. As 
well known from QCD, in view of the property of asymptotic freedom the one-gluon exchange potential 
gives the main contribution at short distances. With the increase of the distance the long-range confining 
interaction becomes dominant. At present the form of this interaction cannot be established in the 
framework of QCD. The most general kernel of qq interaction, corresponding to the requirements of 
Lorentz invariance and of P and T invariance, contains Jl7| ], pSfl scalar, pseudoscalar, vector, axial- vector 
and tensor parts. The analysis carried out in ||, |17j has shown that the leading contributions to the 
confining part of the potential should have a vector and scalar structure. On the basis of these arguments 
we have assumed that the effective interaction is the sum of the one-gluon exchange term and the mixture 
of long-range vector with scalar potentials. We have also assumed that at large distances quarks acquire 
universal nonperturbative anomalous chromomagnetic moments and thus the vector long-range potential 
contains the Pauli interaction. The quasipotential is defined by H 



V(p,q;M) = ui(p)M~P) (ja s iV(k)7f72 + V c v onf (k)T^ + V c b onf (k) j Ul (q)u 2 (-q), (5) 

where a s is the QCD coupling constant, D^ v is the gluon propagator; 7^ and u(p) are the Dirac matrices 
and spinors; k = p — q; is the effective vector vertex at large distances, 



^ (k)=7 ^ + 2m^ 



(6) 
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K is the anomalous chromomagnetic quark moment. 

The complete expression for the quasipotential obtained from (g), (^J) with the account of the rel- 
ativistic corrections of order v 2 /c 2 can be found in Q. In the nonrelativistic limit vector and scalar 
confining potentials reduce to 

KL/M = (1 - e)(Ar + B), V c s onf (r) = e(Ar + B), (7) 

reproducing V™"/ ei (r) = Vf onf + V^ nf = Ar + B, where e is the mixing coefficient. 

All the parameters of our model: quark masses, parameters of linear confining potential A and B, 
mixing coefficient e and anomalous chromomagnetic quark moment k were originally fixed from the anal- 
ysis of meson masses J9| and radiative decays jH^. Quark masses: m\,= 4.88 GeV; m c — 1.55 GeV; 
m s = 0.50 GeV; m U: d= 0.33 GeV and parameters of the linear potential: A= 0.18 GeV 2 ; B= —0.30 GeV 
have standard values for quark models. The value of mixing coefficient of vector and scalar confining 
potentials e = —0.9 has been primarily chosen from the consideration of meson radiative decays, which 
are rather sensitive to the Lorentz- structure of the confining potential JTo| | . Universal anomalous chro- 
momagnetic moment of quarks k = — 1 has been fixed from the analysis of the fine splitting of heavy 
quarkonia 3 Pj states ||. 

Recently, in the framework of RQM, the 1/mg expansion of the matrix elements of weak currents 
between pseudoscalar and vector heavy meson states has been performed Q . It has been found that the 
particular structure of 1/ ttiq corrections up to the second order predicted by HQET can be reproduced 
in RQM only with some specific values of k and e. The analysis of the first order corrections || allowed 
to fix k = —1, while from the consideration of the second order corrections it has been obtained that 
mixing parameter e should be e = —1. Thus HQET, and hence QCD, imposes strong constraints on the 
parameters of the long-range confining potential. The obtained value of e is very close to the previous one, 
determined phenomenologically from radiative decays ifTof and the value of k coincides with the result, 
obtained from the mass spectra ||. Therefore, there is an important QCD and heavy quark symmetry 
motivation for the choice of the main parameters of our model: e = — 1, « = — 1. 



3 Exclusive Semileptonic Decay 
3.1 Form Factors and Decay Widths 

For semileptonic decay B — > A{A*)lvi of pseudoscalar meson B into pseudoscalar (vector) meson A(A*) 
the differential width can be written as 

dY{B -> A(A*)lv;) = -L- | A(B A{A*)l n ) | 2 d$, (8) 

where 

= (2n fS^ (p B -vi- Pvi - PA) J! jSfw. > i = Al,*l, (9) 

Pb is four-momentum of initial meson, p^ is four-momentum of final meson, pi and p Vl are four-momenta 
of lepton and neutrino respectively. 

The relevant transition amplitude looks like 

A(B - A{A*)lm) = (A(A*)lvi\H eff \B) = ^V ab L^, (10) 



where 



V2 

V a b is the CKM matrix element connected with b — > a transition 



Gp 

Beff — ~~/^JhadronJ le P ton ,Hi (H) 
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The leptonic and hadronic currents are denned by 

^=/7„(l- 75 K (12) 

H^ = (A(A*)\a^(l~ l5 )b\B), (13) 

the initial meson B has the quark structure (bq) and the final meson A(A*) has the quark structure (aq). 
The matrix element of hadron current can be expressed in terms of Lorentz-invariant form factors 

a) For CP — > CP transition B — > 

(A(pa)|J^|b(pb)> = U(q 2 )(PA+ P B) f , + f-(q 2 )(PB~PA),; (14) 

b) For CT — > 1" transition 5 — > 



Wa 2 



W) 

4s(<r) 



(A*(p A ,e)|J;?|B(pi,)) = (M A + M B )A 1 (g 2 )e* - — (e*p B )(p A +p B ) A 



M A + M B 



(e*pB)(pB -Pa) v.] (16) 



where g = — Pa; ■Tjf = {a-l^b) and J/f = (07^75^) are vector and axial parts of the weak current; e M 
is the polarization vector of A* meson. 

Since q = pi+p Vl , the terms proportional to q v , i.e. /_ and A3 give contributions proportional to the 
lepton masses and do not influence significantly the transition amplitude, except for the case of heavy r 
lepton, and thus will not be considered. 

The matrix element of the local current J between bound states in the quasipotential method has the 
form @ 

(A\J,(0)\B) = J ^^§4p)r M (p,q)vMq), (17) 

where T M (p, q) is the two-particle vertex function and 'J a,b are meson wave functions projected onto the 
positive energy quark states. 

In the case of semileptonic decays J p = Jhadron;^ — JX — J a ^ s the weak quark current and in order 
to calculate its matrix element between meson states it is necessary to consider the contributions to T 
from Figs. 1 and 2. The vertex functions obtained from these diagrams look like 



rW(p,q) =u a ( Pl ) 1 ^l- l5 )u b (qi)(27r) 3 6(p 2 -q 2 ), (18) 



and 



itVq) = u a (jp x )u q {p 2 ) 7im(1-7i 5 ) »_/ r Tfyfo-qa) 
^ L £b(k)+Sb[Pi) 

q2) i^yT^) 7i ° 7l ^ (1 " 7i5) l u *(ri u M> ^ 

wherek = p!-A; k' = q ± + A; A = p B p A ; e(p) = (m 2 + p 2 ) 1/2 ; 

A (_ } = e{p) - (to 7 ° + 7 °7p) 
2e(p) 

As one can see, the form of relativistic corrections resulting from Tj, (p, q) is explicitly dependent on the 
Lorentz-structure of quark-antiquark potential. 
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Our previous analysis of the semileptonic B — > D{D*) and D — > K(K*) transitions juj was based 
on the assumption that we could expand ( [To| ) up to the order p 2 /m 2 with respect to both b and a 
quarks. This assumption proved to be quite adequate in the case of B — » Dlvi where both b and c quarks 
are heavy. However, the final s quark is not heavy enough. It would be more accurate not to expand 
r^ 2) (p,q) at all, but one should do it in order to perform one of the integrations in fll7| ) . Our statement 
is that more reliable results for semileptonic D — > K decays can be obtained by using p^/e 2 (p) expansion 
instead of p 2 /m 2 in (19). 

It is also necessary to take into account that the wave function of final A meson ^A,p A (p) is connected 
with one in A rest frame ^a.o(p) as follows 

*a, pa (p) = ^yw P j<Wpj*Ao(p), (20) 

where D 1 / 2 (R) is well-known rotation matrix and R w is the Wigner rotation. 

The meson functions in the rest frame have been calculated [^0| by numerical solution of the quasipo- 
tential equation (^). However, it is more convinient to use analytical expressions for meson wave functions. 
The examination of numerical results for the ground state wave functions of mesons containing at least 
one light quark has shown that they can be well approximated by the Gaussian functions 

i*\ 3 " i p 2 



*"<P> = *M,„(p) = ^rj -P^J. m 

with the deviation less than 5%. 
The parameters are 

Pb = 0-41 GeV; (3 K = K , = 0.33 GeV; /3 = 0.36 GeV; 

[3 D = 0.38 GeV; (3 Ds = 0.44 GeV (22) 

In the B meson rest frame equations (|l4|)-(|l6|) can be written in the three-dimensional form 

a) For 0~ — ► 0~ transition B — > Alvi 

(A(p A )\4\B(p B )) = f + (q 2 )(M B + E A ) + f_(q 2 )(M B -E A ), (23) 
(A(p A )\3 v \B( PB )) = (/_(g 2 )-/ + ( 9 2 ))A, (24) 

b) For 0~ — > 1~ transition B — > A*lvi 

V(a 2 ) 

(A*(p A ,e)\jV\B( PB )) = i K l> e 0up<r e*"(pA+PBy(pB-p A r = 0, (25) 

M A + Mb 

(A*(p A ,e)\3 v \B( PB )) = i j^rf V(q 2 )[e*A], (26) 

Ma + Mb 

(A*(p A ,e)\J A \B(p B )) = eg (A 1 (a 2 )(M A + M B ) - A 2 (q 2 ) M ^ Mg (Mb + E A ) 

+MQ 2 ) T7^rT7 (M B -E A )), (27) 

i A 



M A + M B 

(A*(p A ,e)\J A \B(p B )) = A 1 (q 2 )(M A + M B )e Agg m ^ Mb (Mq 2 ) + Ml 2 )), (28) 



where e M = (0, e) is the polarization vector of A* meson in its rest frame, e M is the vector obtained from 
e M by the Lorentz transformation Lj\ 

e M = L^e^. (29) 

The components of e M look like 

eA . A(Ae) . A 

e °=M^' e = e+ M A (E A + M A ) =e + eo E A -TM A -- (30) 
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Equations fl23|), ( P4[) and (|25|)— (p8|) determine form factors /+, /_ and V, A±, A%, A 3 respectively. 

Substituting the vertex functions (|l^) and (|l^) , with the account of wave function transformation ( ^o| ) 
and quasipotential equation (|l]) in the matrix element ( jl7j ) and using eqs. (^3|), ( |24| ) and (p5|)-(p8|) we 



a 2 



{M E 



Imax 

f 

i ll! (IX 



M A ) 

f( 2 )/ 



point 



(1- 

^1 (?roa:c) + £ ^is(9maa;) + (1 
) = ^2 'feLi) + £ ^2s(9maa;) + (1 

) = VW(ql ax )+eVP(qi ax ) + (l 



2 ) 

max J i 



-s)A§(q, 



2 ) 

max J i 

2 ) 

max J ' 



(31) 
(32) 
(33) 
(34) 



get the following expressions at q 2 

/+(<?; 

■^■2 (,Qrnax 

V(q 2 
v x^max 

where /| 2 ] y , Af SiV , A^ v , and are given in Appendix A. In ©-© indeces 

(1) and (2) correspond to the diagrams in Figs. 1 and 2, S and V correspond to the scalar and vector 
potentials of quark-interaction. 

Now our concern is to find q 2 dependence of the form factors. Components of axial and vector currents 
can be expressed in terms of two functions i*i(A) and F 2 (A), A = pb — p A , 

4(A) = F 2 (A), 

A + i[e*A] 



J- (A) 



3 A (A) 



2m a 
(e*A 



L ^i(A), 



(35) 
(36) 



*i(A), 



2m a 
e*F 2 (A). 

Functions F\ and F 2 arise from the lower and the upper components of Dirac spinors 



and are equal to 



Near q 



2 _ 2 

J fit ■ 



Fi(A) = 

F 2 (A) = 
it can be written as 
Fi(A) = 

F 2 (A) = 



2£a(p) 



2m„ 



1/2 



1 



X 



e a (q + A) 4 

e a (q + A) 



£q(q + A) 
2e a (q + . 

1/2 



1/2 



Ea 
M a 



2e a (q+ A) 
V2(U 



Ea 
M a 



A 2 /M 2 ) l/2 



(l + A 2 / m 2 . 




y/l + A 2 /m 2 a 

\ 1/2 
1 \ 



1/2 • 



+ A 2 /m 2 a 



1/2 



(37) 
(38) 

(39) 

(40) 
(41) 

(42) 
(43) 



The dependence of the formfactors on the momentum transfer is fixed by extrapolating their behavior 
(A = 0) point over the kinematically allowed region 



near q 2 = q 2 max 



f+(A) = /+(0)/(A)F 1 (A), 

At(A) = Ax(0)7(A)F 2 (A), 

A 2 (A) = A 2 (0)I(A)F 1 (A), 

V(A) = V(0)I(A)F 1 (A), 



(44) 
(45) 
(46) 
(47) 



(i 



where 



Introducing the variable 



M 2 +M A -q 2 

vavb = B n , T * ( 49 ) 



2M A M B 

where v A and v B are meson velocities, and taking into account that 

9 , „ (Ml + Ml - q 2 ) 2 , , 

A 2 = ( PB - Pa) 2 = 1 ai\/t2 -M 2 =M 2 ( W 2 -1), (50) 



4M1 



we can rewrite (|44|)— (47) in the form 



1/2 



/ + H = /+(!)/(«,)(_! , (51) 
.!,(«•) = .l l , 1 ,A(„, ) (^V /a , (52) 



1/2 



A 2 H = A a (l)/(«;) — — , (53) 
\u> + 1/ 

/ 2 \ 1/2 

n») = V(l)I(w) — — , (54) 
\u> + 1/ 

Substitution of the Gaussian wave functions (^l]) in (E§) results in 

' w =^(-?i^k^> (1) - (55) 

where A =< e 9 > is a mean value of light quark energy inside meson. In our model A corresponds to 
HQET parameter A = M — toq, which determines the energy carried by light degrees of freedom, and is 
found to be g: A = 0.54 GeV. 

In the limit of infinitely heavy b and a quarks the w dependence of eqs. (|5l|)-(|55|) is determined by 
the Isgur-Wise function of our model 

1/2 ' K 2 w-\ 



« w > = U+iJ exp \~W^+i) (56) 

and the ratios of form factors satisfy all constraints imposed by HQET . 

Using (H)-(H) and (??)-(??) we have calculated form factors for B -> D{D*)lv h D -> K(K*)lv t 
and D s — > iplvi exclusive decays. The results obtained in our model for D — ► K(K*)lvi are compared with 
appropriate experimental data and various model predictions in Tabic 1. New values of form factors for 
D — > K(K*)lvi are somewhat larger then our previous results JTJ] because of more consistent relativistic 
treatment of s quark and slight change in the value of the mixing coefficient e. Note that while the other 
potential models agree with the experimental determination of V(0), but fail to predict Ai(0) and ^2(0), 
our model predicts correct values of Ai(0) and ^b(O), but gives too low value of V(0). The reason for that 
is not clear. The contribution from form factor V(0) in the total width is kinematicaly suppresed. So, 
despite the above mentioned descripancy, we have got the D — ► K*lv\ width in accord with experimental 
data. 

The ratios of form factors R 2 = A 2 (0)/A 1 {0) and R v = V(0)/A x (0) arc given in Tabic 2. 
The obtained branching ratios are 

B(D° -> K*~e+v e ) = 1.9%, for t d o = 0.415 x 1Q- I2 s; 
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B{D+ -» K*°e+v e ) = 4.9%, for T£)+ = 1.060 x 10- 12 s; 
to be compared with the experimental average data [^2f 

B ex P ( D o _> K *- e + v ^ = (2.0 ± 0.4)%, 

B ex P(D+ -> K*°e+v e ) = (4.8 ± 0.4)%. 
The ratio — * K*ev e )/T{D — > Kev e ) and the ratio of the longitudinal and transverse decay widths 
Tl/I^t also agree well with experiment (see Table 3). 

For D s — > iflvi decay form factors our predictions are 

Ai(0) = 0.63, A 2 (0) = 0.35 and V(0) = 1.06. 
As experiment provides us with R 2 and Ry ratios for D s — > tplvi, we compare it with predicted values 
in Table 4. It is not clear why experimental ratio R 2 for D s — -> iplvi so differs from that of D — > K*lvi. 
In RQM we get approximately equal ratios R2 for both decays, because general structure and the signs 
of the potential-dependent corrections in (J3l"|)~(|34l) are almost the same. It can be expected that the 
experimental results for D s — > iplvi form factor ratios will change in the future. Anyway, the experimental 
uncertainties are still too large to conclude that there is a serious descripancy between RQM and the 
experimental data in this case. 

For D s — ► tplvi branching ratio we have B(D S — > <plvi) — 2.5%, while experiment gives B exp (D s — > 
tplvi) = (1.88 ±0.29)% [||. 

In B — > D*lvi decay, since both b and c quarks are heavy, relativistic corrections are not so significant, 
but the Lorentz- structure of quark-antiquark potential has an important influence on the values of form 
factors. We have found our results for R2 and Ry to be in a good agreement with the experimental 
data [^3) and HQET-based predictions [^5). Measurments and predictions for the ratios of the form 
factors for B — > D*lvi, evaluated at q 2 = q 2 lax , are shown in Table 5. We have obtained the following 
results for B — > D*lvi and B — > Dlvi branching ratios 

B(B -> D*lvi) = 33.8 x \V bc \ 2 , B(B -> Din) = 19.8 x \V hc \ 2 , for t b o = 1.5 x 10- 12 s. 
It should be compared to the experimental data 



B(B° -> D-e+u e ) = (2.0 ± 0.7 ± 0.6)% ARGUS 26 



B(B° D-e+u e ) = (1.8 ± 0.6 ± 0.3)% CLEO- I [27 



B{B° -> L>*-e+iy e ) = (4.7 ± 0.5 ± 0.5)% ARGUS |28 



5(5° -► D*-e+v e ) = (4.0 ± 0.4 ± 0.6)% CLEO- I |7 
As a result we can extract the value of CKM matrix element V c b 

\Vcb\ = 0.036 ±0.004. 



3.2 Differential Distributions 



The differential decay rate \21 can be expressed in terms of two dimensionless variables x = Ei/Mb and 
w = vavb, where E\ is the lepton energy 



d 2 T 

dxdw 



\Vg b \ 2 G 2 F M% 
32tt 3 



(2M A 
V M B 



+2W 2 {w)( K p + (1 - 2xf - -j^-wQ. - 2xj) 



1 - pj (W^w) - 2W 3 (w)(l -2x- j±w) 
2 Ma 



(57) 



Here W\ 2,3 ( w ) are connected with semileptonic form factors 
a) For 0~ — * _ transition 



Wi(w) 
W 2 {w) 



0, 

2M A 
Mb" 



I/+HI 5 



W 3 (w) = 0, 



(58) 
(59) 
(60) 



s 



b) For — > 1 transition 



W = 2M- a \ + M~b) A M ~ M A + M B V {W) { 1 + P ~ U^ W ) 

+2Ai(w)A 2 (w) - «) + -g^^)^ - 1), (62) 

w 3 H = -i/iiW^H. (63) 

M B 

The kinematiclly allowed region is presented in Fig. 3 where the lower bound curve w m {x) has the 
following shape 

, x m b M A 1 . , 

Wm (x ) = — 2_(i_2a;) + — - . (64 

y ' 2M A ' 2Mb l-2x 1 1 



The analytical expression for dT/dx distribution depends on q 2 behavior of form factors. Using (|5l|)-( 
we obtain after the integration over w 



dT _ G 2 F \V cb \ 2 M% 



dx 32tt 3 
where 



T^K x {x) +sinh{a X -{x))e- ax+{x) K 2 (x) + K 3 (x) / dwe~ a ^ 

J w m (x) 



(65) 



x {X ) = AMaM b < l - R ) (66) 

(Ma + Mb^i- 2M » , xa-R) 



\M B + Ma J 1 2M b T -1 _ m 

p _ M\ 4A 2 _ Aff + M\ 

1-2*' p -M 2 > a ~ P 2 Ma +(3 2 Mb > Wo ~ 2M A M B - 

Functions Ki^ t 3(x) take different forms for 0~ — > 0~ and 0~ — » 1 — decays and are given in Appedix B. 

The (l/T)(dT/dx) distributions for B — ► D{D*)lvi, D — ► K{K*)lvi and D s — > sfl/z/; decays are shown 
in Figs. 4- 6. All curves are normalised by the corresponding decay width T, i.e. the area under each 
curve is equal to 1. 



4 Conclusion 

Using the quasipotential approach, we have obtained the expressions for semileptonic decay form factors 
with the consistent account of relativistic effects. This account includes more careful consideration of 
the s quark contribution than it was in our previous work |l4| and results in the small shift in the 
values of D — > K{K*)lvi decay form factors, which are in a good agreement with measurments. Our 
model provides more accurate values for the ratio A2(0)/yli(0) in the D — > K{K*)lvi decay and for the 
ratio T(D — > K*lvi) /T{D — > Klvi) in comparison with the other models. We have also calculated form 
factors and branching ratios for B — ► D(D*)lvi, D — > K(K*)lvi and D s — ► cplvi. The extracted value of 
\Vcb\ = 0.036 ± 0.003 is lower than the previous one Q because of the changes in the form factors as 
well as in the B meson lifetime. We should emphasize that in order to get reliable results for D meson 
semileptonic decays it is necessary to take into consideration all possible relativistic effects, including the 
transformation of meson wave function ( |20| ) from the rest frame, which is ignored in many quark models. 
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The proposed q 2 dependence of the form factors is used for the determination of differential semilep- 
tonic distributions in the case of pseudoscalar and vector final states. 

It should be noted that obtained expressions for semileptonic form factors are valid for all B and 
D meson decays, except the decays into mesons containing two light quarks (n, p mesons), where one 
cannot expand in either p 2 /m 2 or p 2 /e 2 at q 2 = point in the vertex function Jl9|). The solution of 
this problem is proposed in [Q. 

The analysis has shown that Lorentz-structure of quark-antiquark potential plays an important role 
in heavy meson semileptonic decays. We have got experimentally motivated and HQET based arguments 
to conclude that the confining potential has predominantly Lorentz-vector (with Pauli-term) structure 
e = — 1. Assuming also long-range anomalous chromomagnetic moment of quark to be n — — 1 we have 
obtained satisfactory description of all considered B and D semileptonic decays. 

We argue that small number of parameters, most of which were fixed previously, and an agreement 
with HQET for the structure of leading, subleading and second order terms in 1/mg expansion make 
RQM a reliable tool for the investigation of heavy meson physics. 

In this paper we have restricted ourselves with 0~ — > 0~ and 0~ — > 1 _ transitions. One more 
practically important case of the decay into P wave final state (i.e. B — > D**lv{) will be considered in 
the future. 
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Appendix A. Exclusive Semileptonic Decay Form Factors at q 2 = 

Qmax P° int 

arXiv note - added Dec99 - TeX source for this appendix was corrupt on initial submis- 
sion, changed to verbatim environment so that the rest of the paper can be processed. 

\medskip 

\begin{eqnarray} 

\label{35} 

\nonumber 

f_+"{(l)X\qm)&=&\sqrt{\frac{M_AHM_B}}\int\frac{d-3{\bf p}}{ (2\pi) ~3} 
\bar\Psi_A({\bf 

p})\left (\frac{\varepsilon_a+m_a}{2\varepsilon_a}\riglit) "{1/2} 

\Bigg[l+\frac{M_B-M_A}{\varepsilon_a+m_a} 

-\frac{{\bf p}~2}{8}\bigg(\frac{l}{m_b~2} 

-\f rac{4}{m_b (\varepsilon_a+m_a) }\\ 

\nonumber 

&&+\frac{M_B- 

M_A}{\varepsilon_a+m_a}\Big(\f rac{l}{m_b~2}+\f rac{4}{\varepsilon_a 
(\varepsilon_a+m_a)}+\frac{2}{3\varepsilon_am_a}\Big)+\f rac{4}{3}\f rac 
\nonumber 

&&\times\Big(\irac{l}{\varepsilon_a+m_a}- 
\frac{l}{2m_b}\Big)\f rac{\varepsilon_a 
-\varepsilon_qM\varepsilon_a \varepsilon_q} \bigg) 
+\frac{M_B-M_A}{3}\bigg(\frac{l>{\varepsilon_a+m_a}+ 
\frac{lH2m_b}\\ 

&&-\frac{{\bf p}~2}{8m_b~2}\Big( 

\frac{l}{\varepsilon_a+m_a}+\frac{3}{2m_b} 

\Big) \Bigg) \frac{\varepsilon_q}{M_A} ({\bf p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p»,W 

\label{36} 

\nonumber 

f_{+S}-{(2)}(\qm)&=&\sqrt{\frac{M_A}{M_B}}\int\frac{d-3{\bf 

p>K(2\pi)-3> 

\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a}\righ.t) "{1/2} 
\Bigg[-\f rac{M_B-M_A}{2\varepsilon_a}\f rac{M_B-\varepsilon_b- 
\varepsilon_q} 
{\varepsilon_a}+ 

\frac{M_B-M_A}{12\varepsilon_a}{\bf p}~2\\ 
\nonumber 

&&\times\Big(\frac{l}{\varepsilon_a"2}+\frac{l}{m_b~2}\Big)- 

\frac{M_B-M_A}{12}\Big(\frac{l}{\varepsilon_a-2}+\frac{l}{m_b-2}\Big) 

\big(M_B+M_A-\varepsilon_b-\varepsilon_a-2\varepsilon_q\big) \\ 

&&\times\frac{\varepsilon_q}{M_A} ({\bf p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p»,W 

\label{37> 

\nonumber 

f_{+V}"{(2)}(\qm)&=&\sqrt{\frac{M_A}{M_B}}\int\frac{d-3{\bf 

p}}{(2\pi)~3} 

\bar\Psi_A({\bf 

p})\left (\frac{\varepsilon_a+m_a}{2\varepsilon_a}\riglit) "{1/2} 
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\Bigg[\frac{M_B-M_AM\varepsilon_a}\f rac{{\bf 
p}~2Ml2}\bigg( (l+\kappa) 

\Big(\frac{l}{\varepsilon_a~2}Afrac{l}{m_b~2}\Big)\\ 

\nonumber 

&&- 

\f rac{l}{\varepsilon_q}\Big(\f rac{2}{\varepsilon_a+m_a}+\f rac{l}{m_b} 
\Big)\bigg)+ 

\frac{M_B-M_A}{12}(l+\kappa)\Big(\frac{l}{\varepsilon_a-2>- 
\frac-ClMm_b~2} 

\Big)\big(M_B-M_A-\varepsilon_q+\varepsilon_a\big)\\ 
\nonumber 

&&\times\f rac{\varepsilon_q}{M_A} ({\bf p}\overlef tarrow 
{\irac{\partialH\partial{\bf p}») 

+\f rac{M_B-M_A}{6\varepsilon_q}\Big(\f rac{l}{\varepsilon_a+m_a} 
+\f rac{l}{2m_b}\Big) \big(M_B+M_A-\varepsilon_b-\varepsilon_a- 
2\varepsilon_q\big) \\ 

&&\times\f rac{\varepsilon_q}{M_A} ({\bf p}\overlef tarrow 
{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p»,W 
\nonumber 

A_l-{(l)}(\qm)&=&(M_A+M_B)\sqrt{4M_AM_B}\int\frac{d-3{\bf 

p>K(2\pi)-3} 

\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a}\right) "{1/2} 

\bigg(l-\frac{{\bf p}-2>{8>\Big(\frac{l}{m_b"2}\\ 

&&+\f rac{4}{3 (\varepsilon_a+m_a)m_b>\Big) \bigg) \Psi_B ({\bf p» , \\ 

\label{39> 

A_{lS}-{(2)}(\qm)&=&A_{lV}-{(2)}(\qm)=0,\\ 

\label{40> 

\nonumber 

A_2-{(l)}(\qm)&=&\frac{l}{(M_A+M_B)\sqrt{4M_AM_B»\int\frac{d-3{\bf 
P» 

{(2\pi)~3}\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a> 

\right)-{l/2}\Bigg[\bigg(l+\frac{M_A}{M_B}\bigg)\\ 

\nonumber 

&&\times\bigg(l 

-\frac{{\bf 

p}"2}{2}\Big(\frac{l>{4m_b~2}+\f rac{l}{3m_b(\varepsilon_a+m_a)> 
\Big)\bigg)- 

\frac{2M_A-2}{M_B(\varepsilon_a+m_a)}\bigg(l-\frac{{\bf p}~2}{8} 

\Big(\frac{l}{m_b~2}\\ 

\nonumber 

&&+\f rac{4}{\varepsilon_a(\varepsilon_a+m_a) }+ 

\f rac{2}{3\varepsilon_am_b}\Big) 

-\frac{{\bf p}~2H6M_A\varepsilon_a} 

\Big(l+\f rac{\varepsilon_a+m_a}{2m_b}\Big)\bigg)\\ 

\nonumber 

&&-\f rac{2M_A\varepsilon_q}{(\varepsilon_a+m_a)M_B}\bigg( 
\frac{l}{3}\Big(l+\frac{\varepsilon_a+m_a}{2m_b> 
-\frac{{\bf p}~2H8m_b~2} 

\Big(l+\frac{3(\varepsilon_a+m_a)M2m_b}\Big)\Big)\\ 
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&&+\frac{M_B-M_AH3m_b}\Big(l-\frac{3{\bf p}~2}{8m_b~2}\Big) \bigg) 
({\bf p}\overlef tarrow 

{\frac-C\partialH\partial-C\bf p»»\Bigg] \Psi_B({\bf p»,W 

\label{41> 

\nonumber 

A_{2S}~ { (2) } (\qm) &=&\f rac{l}{ (M_A+M_B) \sqrt{4M_AM_B}}\int\frac{d~ 3{\bf 
{(2\pi)'3}\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a} 

\right) ~{l/2}\f rac{M_A"2}{\varepsilon_aM_B}\Bigg [-\f rac{M_B- 

\varepsilon_b 

-\varepsilon_q}{\varepsilon_a}\\ 

\nonumber 

&&+\frac{{\bf 

p}"2}{6}\Big(\frac{l}{\varepsilon_a-2}+\frac{l}{m_b-2}\Big)- 

\f rac{\varepsilon_a}{6}\big(M_B+M_A-\varepsilon_b-\varepsilon_a- 

2\varepsilon_q\big) 

\Big(\f rac{l}{\varepsilon_a~2}+\f rac{l}{m_b~2}\Big) \f rac{\varepsilon_q 
{M_A}({\bf p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p»,W 
\\ 

\label{42> 
\nonumber 

A_{2V}~ { (2) } (\qm) &=&\f rac{l}{ (M_A+M_B)\sqrt{4M_AM_B}}\int\f rac{d~3{\bf 
{(2\pi)-3}\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a} 

\right)'{l/2}\Bigg [\f rac{MJT 2H\varepsilon_aM_B}\bigg(\f rac{{\bf 

p}~2}{6> 

(l+\kappa)\\ 

\nonumber 

&&\times\Big(\frac{l}{\varepsilon_a"2}-\frac{l}{m_b~2}\Big)+ 

\f rac{\varepsilon_a}{6}\big(M_B-M_A-\varepsilon_b+\varepsilon_a\big) 

(l+\kappa)\Big(\frac{l}{\varepsilon_a~2}-\frac{l>{m_b"2}\Big) 

\f rac{\varepsilon_q}{M_A} ({\bf p}\overlef tarrow 

{\f rac{\partial>{\partial{\bf p»» \bigg) \\ 

&&-\irac{{\bf p}~2}{6\varepsilon_q} 

\nonumber 

\Big(\f rac{2}{\varepsilon_a+m_a} 
+\frac{l}{m_b}\Big)+ 

\f rac{\varepsilon_a}{6M_A}\Big(\f rac{2}{\varepsilon_a+m_a}+\f rac{l}{m_ 
&&\times\big(M_B+M_A-\varepsilon_b-\varepsilon_a-2\varepsilon_q\big) 
({\bf p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p»,W 

\label{43> 

\nonumber 

V-{(l)}(q-2_{max»&=&\frac{l}{M_A+M_B}\sqrt{\frac{M_A}{M_B»\int\frac{ 
-C(2\pi)~3}\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a} 

\right)"{l/2}\frac{l}{\varepsilon_a+m_a}\Bigg[l-\frac{{\bf p}~2}{8} 

\Big(\frac{l}{m_b~2}\\ 

\nonumber 

&&+\f rac{4}{\varepsilon_a(m_a+\varepsilon_a)> 
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-\f rac{2}{3\varepsilon_am_b}\Big)- 

\irac{{\bf p}~2}{12M_A>\Big(l+\frac{\varepsilon_a+m_a}{2m_b>\Big) 
\f rac{\varepsilon_a+\varepsilon_qM\varepsilon_a\varepsilon_q}- 
\frac{{\bf p}~2}{12M_A\varepsilon_a}\Big(l- 
\frac{\varepsilon_a+m_a}{2m_b}\Big) \\ 

&&+\frac{l}{3}\Big(l-\frac{\varepsilon_a+m_a}{2m_b}-\frac{{\bf 
p}~2M8m_b~2} 

\Big(l+\frac{3(\varepsilon_a+m_a)>{2m_b>\Big)\Big) 

\f rac{\varepsilon_q}{M_A> ({\bf p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p»,W 

\label{44> 

\nonumber 

V-{(2)}_S(q-2_{max})&=&\frac{l}{M_A+M_B}\sqrt{\frac{M_A}{M_B}}\int\fra 
{(2\pi)-3}\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a> 

\right) "{l/2}\f rac{l}{2\varepsilon_a}\Bigg[-\f rac{M_B-\varepsilon_b 

-\varepsilon_q}{\varepsilon_a}\\ 

\nonumber 

&&+\frac{{\bf p}~2M6}\Big(\frac{lM\varepsilon_a-2}- 
\frac{lMm_b~2}\Big)- 

\f rac{\varepsilon_a}{6}\big(M_B+M_A-\varepsilon_b-\varepsilon_a- 
2\varepsilon_q\big) \\ 

&&\times\Big(\frac{l}{\varepsilon_a~2}— 
\f rac{l}{m_b"2}\Big)\f rac{\varepsilon_q> 
{M_A}({\bf p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p}),\\ 

\label{45> 

\nonumber 

V-{(2)}_V(q-2_{max»&=&\frac{l}{M_A+M_B}\sqrt{\frac{M_A}-[M_B»\int\fra 
{(2\pi)-3}\bar\Psi_A({\bf 

p})\lef t (\f rac{\varepsilon_a+m_a}{2\varepsilon_a} 

\right)-{l/2}\frac{l}{2\varepsilon_a}\Bigg[\Big(\frac{{\bf p}"2}{6> 

+\f rac{\varepsilon_a}{6}\\ 

\nonumber 

&&\times\big(M_B-M_A-\varepsilon_b+\varepsilon_a\big) 
\times\f rac{\varepsilon_q}{M_A}({\bf 
p}\overlef tarrow{\f rac{\partial}{\partial 
{\bf p»»\Big) 

(l+\kappa)\Big(\frac{l}{\varepsilon_a~2}+\frac-[l}{m_b"2}- 

\f rac{l}{\varepsilon_q}\Big(\f rac{2}{\varepsilon_a+m_a}+ 

\frac{l}{m_b>\Big)\Big)\\ 

\nonumber 

&&-\frac{{\bf 

p}~2\varepsilon_a}{6\varepsilon_q}\Big(\frac{2}{\varepsilon_a+m_a} 

-\frac{l}{m_b}\Big) 

+\f rac{\varepsilon_a}{6M_A} 

\big(M_B+M_A-\varepsilon_b-\varepsilon_a-2\varepsilon_q\big) \\ 
&&\times\Big(\frac{2}{\varepsilon_a+m_a}-\frac{l}{m_b>\Big) ({\bf 
p}\overlef tarrow 

{\frac{\partial}{\partial{\bf p»»\Bigg] \Psi_B({\bf p» , 
\end{eqnarray} 
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here $({\bf p}\overlef tarrow{\partial/\partial{\bf p}})$ acts 
to the left on the wave function $\bar\Psi_A({\bf p})$. 
In the limit ${\bf p}~2/m~2\to 0$ the above form factors reduce to the 
standard expressions, obtained in the nonrelativistic quark models. 



Appendix B. Functions K\ ,2,3 (#) for [l/Y)dY / dx Differential Distri- 
butions 

a) In the case of - — > 0~ transition 

Ki(x) = Hgl a ,(l-2 I )(1-R)|/ + (1)| 2 , (68) 

" ^7( 1 + ^) 2<1 -«l /+(1 »l 2 ' (69) 



* 3 (aO = J ^[p+{l-2xY + 1 ^{l-2x)- 1 ^a{l-2x))\U{l)\\ (70) 
b) In the case of - — > 1~ transition 



ifi(x) = ^x(l-i?)(G 1 (a;)+ a G 2 ( a; ) + ^ 2 G 3 ( a; )) + i(G 2 (x) + ^aG 3 ( a; )) 

/ , w M Bn 2 Mo 9/ , 9 \ 1 , X /3M B , 
x (,(1 - r)(l + jf ) - jfi/d - *) 2 ) + 3^(^(1 - R) 



(M A + M B ) 4 



A 1V± B 



3_Ml ( M A + M B y ( M A + M B y x 

+ Mi V 2 Ma Mb ) [ ' V 8MiM| ) 1 j J' 1 j 

^ (a;) = - (M m^ b)2 + aG2(a;) + r 2G3(x)) 

(Ma + Mb)* 2 (Ma + Mg) 6 
MP^( G2(X) + r° 3 ^ 12M1M3 (72) 

£T 3 (z) = -G 4 (.T)-2aGi(x)-2a 2 G 2 (a;)-|a 3 G3(x), (73) 
where Gi, 2,3,4^) depend on the values of form factors at w = 1 point 



(M^ + Mb) 2 ^ w V Mb/ mb 

2 



MB A,,l M2 (l) l l + -^ jl l-2 I+ -^j , (74) 
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aw . _^^ l(1M2(1)(1 _ 2l) _l^ Al(1)ni)+ _i^_(4 (1)(1 _ 2l) 

-^)) + W( 1+ K)"^>(- 1+ ^- a "). 



G4(x) = (mXT^ f(1 H 1 + m^J Vi - 2 - i+ m^J 
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Table 1. Theoretical predictions and experimental data for form factors in D — > Klvi and D — > D*lvi. 



Ref. 


/+(0) 


Vb(0) 


Ai(0) 


A 2 (0) 


Exp. Average |22j 


0.75 ±0.02 ±0.02 


1.1 ± 0.2 


0.56 ±0.04 


0.40 ±0.08 



Theory 

RQM 0.73 0.62 0.63 0.43 

ISGW H 0.82 1.1 0.8 0.8 

BSW [0] 0.76 1.3 0.88 1.2 

AW 0.7 1.5 0.8 0.6 

BKS||| 0.9 ±0.08 ±0.21 1.4 ±0.5 ±0.5 0.8 ±0.1 ±0.3 0.6 ±0.1 ±0.2 

LMMS H 0.63 ± 0.08 0.9 ±0.1 0.53 ± 0.03 0.2 ± 0.2 

BBD § 0.6 1.1 0.5 0.6 



Table 2. Calculated and measured D 
nO)Mr(O). 



K*lvi form factor ratios i? 2 = A 2 (0)/Ai(0) and i?y = 



Ref. 



H 2 (0) 



i?y(0) 



Experiment 
E691 |0[ 
E653 || 
E687 pl 



0.0 ±0.5 ±0.2 
0.82t°;^±0.11 



2.0 ±0.6 ±0.3 



2.00 



+0.34 
-0.3 2 



±0.16 



0.78 ±0.18 ±0.10 1.74 ±0.27 ±0.28 



Theory 
RQM 
ISGW 
BSW 
AW 
BBD~P 



0.68 
1.0 
1.36 
0.75 
1.2 



0.98 
1.37 
1.48 
1.87 
2.2 



Table 3. The ratios T(D — > K*lv{) /T(D — » Klv{) and T^/IY in comparison with the experimental data. 



Ref. 



Y{K*)/T{K) 



T l /Tt 



RQM 

Experiment 
E691 |0| 
E653 [31 1 
E687 [ p[ 
CLEC^p3| 
CLEO [33 1 



0.65 



0.60 ±0.09 ±0.07 
0.65 ±0.09 ±0.10 



1.05 



^±0.3 
1.18 ±0.18 ±0.08 
1.20 ±0.13 ±0.13 



1. 
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Table 4. Measured and calculated ratios of form factors in D s — > 



Ref. 


i? 2 (0) 


Rv(0) 


Exp. Average |22fl 


1.8 ±0.5 


2.0 ±0.7 


CLEO >jp 


1.4 ±0.5 ±0.3 


0.9 ±0.6 ±0.3 


E653 p6f 


2.1±g;f ±0.2 


2-3^ ±0.4 


Theory 






RQM 


0.55 


0.94 


BKS | 


2.0 ±0.19 ±0.23 


0.78 ±0.08 ±0.15 


LMMS @ 


1.65 ±0.2 


0.33 ±0.36 



Table 5. Predicted and measured ratios of form factors in B — ► D*lvi at q 2 = qf r 



Ref. 






Experiment 






CLEO pi a) 


1.02 ±0.24 


1.07 ±0.57 


CLEO H b) 


0.79 ±0.28 


1.32 ±0.62 


Theory 






RQM 


1.16 


1.74 


ISGW 11 


1.14 


1.27 


WSB |f 


1.06 


1.14 


HQET-based || 


1.26 


1.26 


HQET-based §5 


| 1.14 


1.74 



Figure Captions 

Fig.l Lowest order vertex function. 

Fig. 2 Vertex function with account of the quark interaction. Dashed line corresponds to the effective 
potential (5). Bold line denotes the negative-energy part of quark propagator. 

Fig. 3 Allowed region for semileptonic B — > A(A*)lvi decay in terms of the variables w and x. Lower 
bound curve W m (x) is determined by (64), upper bound is w a — [M\ + M b )/{2MaM b ). 



Fig. 4 (1/T) (dT/dx) for B — > Dlvi and B — > D*lvi. Absolute rates dT/dx can be obtained by using 
T(D) = 1.71 x lO^s- 1 and T(D*) = 2.92 x lO^s" 1 for t b o = 1.5 x 10~ 12 and \V bc \ = 0.036. 

Fig. 5 (l/T)(dT/dx) for D — > Klvi and D — » K*lvi. Absolute rates dT/dx can be obtained by using 
r(iT) = 6.68 x lO 1 ^ 1 and T(K*) = 4.34 x lO 10 ,^ 1 for t d « = 0.415 x 10" 12 s , t d + = 1.06 x 10~ 12 s. 

Fig.6 (l/r)(dT/da;) for D s -> </>Z^. Absolute rate dr/dx can be obtained by using T = 5.42 x lO^s" 1 
for td = 0.47 x 10~ 12 s. 
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